y = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3637, 3507, 3356, 3207, 3008, 2237, 422, 51, 0, 0, 0, 0, 0, 0, ];
x = 50:50:2000;
plot(x, y,'-x')
hold on
y = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3501, 3451, 3401, 3351, 377, 0, 0, 0, 0, 0, 0, ];
plot(x, y, '-*')
hold on
y = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3501, 3446, 443, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y, '-+')
hold on
y = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3550, 921, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ];
plot(x, y, '-^')
hold on
y = [4951, 4901, 4851, 4801, 4751, 4701, 4651, 4601, 4551, 4501, 4451, 4401, 4351, 4301, 4251, 4201, 4151, 4101, 4051, 4001, 3951, 3901, 3851, 3801, 3751, 3701, 3651, 3601, 3551, 3501, 3451, 3401, 3350, 3243, 3107, 2950, 2607, 657, 49, 0, ];
plot(x, y, '--')


